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ABSTRACT 



The acoustic impedapice of several Helmholcz resonators 
were measured using two different techniques. The first 
technique employed the conventional Standing Wave Tube 
method. The second was the Dual-Microphone Transfer Function 
method of J.Y. Chung and D. A. Blasar. The calculability of 

the frequency dependent Helmholtz resonator impedance 

allowed the methods to be compared for both absolute accu- 
racy and relative precision. The average relative precision 
for the SWT technique was 4.2% while the average for the 
dual-microphone technique was 18.4%. 
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3^* INTRO DUCTI3N 



Acoustic iapsdance tubes have bean used for nany years 
to obtain ccaplex acoustic impedances, reflection coeffi- 
cients and phase infcrnacion for many different types of 
materials. The techniques used winn the impedance rubes have 
been many and varied. This study involves the comparison of 
the dual- microphone transfer function technique used by 
Chung and Blaser [Ref. 1] and the standing wave tube (SWT) 
technique to determine their absolute accuracy and relarive 
precision. This work was undertaken as a prelude to the 
extension of computer controlled impedance measuremenr tech- 
niques ro water. 

A. IHPEDANCE HEASDREMENT TECHNIQUES 

Many different impedance measuremenr techniques have 
been used over the years. Beranek [Ref. 2] divides the 
different techniques into three general categories. Surface 
Methods, Transmission Line Methods and Comparison Methods. 
This study looked closely at only two of the techniques, 
however, some of the other methods are worth mentioning as 
background . 

^ • Su r fac e Meth ods 

These methods involve the measurement of the sound 
pressure at a point on the surface of inrerest and the 
particle or volume velocity at the surface. Cook [Ref. 3] 
did this using a short tube of a given radius and length and 
a clamped-edgs diaphram. He was then able to calculate the 
volume velocity and the average pressure over the surface of 
the sample. Beranek [Ref. 4] presented a useful theory for 
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predicring the normal acoustic impsiancs of porous maosrials 
in terras ct their flew resisrance, porosity and sffactive 
dynamic mass of the enclosed gas. 

Trans mission M athod s 

By far, the most popular methods of impedance meas- 
urements fall into this category. The standing wavs tube and 
dual- microphone techniques in this study are included. They 
will be examined in detail later. Another technique involves 
the measurement of the pr essure-fr equency resonance curve 
and determining the curve width and resonant frequency. 
Harris [Hef. 5] [Eef. 6] discusses two of these methods 
which are closely related to each other. 

Finally, one additional method currently in use by 
the O.S.Navy for making acoustical measurements deserves 
mentioning. [Ref. 7] This system is digital and is used for 
underwater acoustical measurements from 1 Hz. to approxi- 
mately 2 MHz. This digital measuring system employes a 
buffered video analog-to-di gital converter and was primarily 
designed to perform acoustical measurements using the Preny 
Method. [Ref. 8] [Ref. 9] The system uses a stable signal 
and sampling source and 12-bit, 5 MHz., analog- to-digital 

converter with a 4K buffer memory. A PDP-11/23 serves as a 
controller as well as a signal processor for averaging and 
amplitude and phase estimation. 

3. Compar ison M e thods 

These methods mainly involve the use of acoustical 
bridges or reliable acoustical standards. Steward first 
introduced [Ref. 10] the simplest and most rapid to use. 
Later refinements followed [Ref. 11] but suitable variable 
standards of acoustical resistance and reactance were diffi- 
cult to obtain. The use of reliable acoustic standards, 
rigid wall and an eighth-wavelength tube, was reported by 
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Flandsrs (Rsf. 12] ir. 1932. Howsv~l, naatheaatical calcula- 
tions vers difficult. 

In uhis study, the transmission line method and the 
comparison method are combined by introducing a new ’’stan- 
dard”, the Helmholtz resonator. This standard has the 
advantage of having an easily calculable reactance which is 
absolute to within one adjustable parameter. This parameter 
is fixed by the Helmholtz resonance frequency which is 
easily measured to a high degree of precision and absolute 
accuracy . 

B. HISTORY OF DUAL- MICROPHONE MEASHRING TECHNIQUES 

Clapp and Firestone [Ref. 13 ] first mentioned the use of 
two microphones in an article published in 19hl. m this 
article it was recognized that their equipment could be used 
to measure the absorption coefficient of materials. Later, 
Shultz [Ref. 14 ] published an article on the acoustic watt- 
meter and its application for measuring absorption coeffi- 
cients. During this time the use of the acoustic wattmeter 
fcr impedance measurements was little used since the 
standing wave tube (SWT) technique was the most fevered 
method. 

In 1977 Seybert and Ross [Ref. 15] published an article 
which discussed an experimental technique using two micro- 
phones and random noise excitation. This technique used 
digital methods and the auto- and cross- spectra to determine 
the normal complex reflection coefficient cf a sample termi- 
nating a tube. Finally, in 1900, Zhung and Blaser [Ref. 1] 
published their article discussing the transfer function 
technique for measuring in-duct acoustic properties. This 
technique also used digital methods and the auto- and 
cross-spectra to determine the normal complex reflection 
coefficient. However, the mathematics did not involve the 
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matrix-inversion used by Seybart and Ross and as a resuir 
the computation time was significantly reduced. 

A two microphone technique utilizing pure tone rather 
than random excitation and analog rather than digital 
processing is described by Blliot. [Ref. 16] 

C. P0RPCS2 OF STUDY 

The above methods have been used to maasure numerous 
materials of differing acoustic properties. However, the 
properties of these materials were not known or calculable 
to any reasonable accuracy. The aim of this study was to 
determine the absolute accuracy of the dual-microphone tech- 
nique and its relative precision compared to the SWT tech- 
nique as a prelude to the development of an impedance 
measurement apparatus for use in water. A Helkholtz reso- 
nator was used to accomplish this comparison. 
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II 



• COaPAR ISON OF MEISOREMENT TECHNIQUES 

Previous studies of impedar.ca measurement techniques 
involved the use of materials for which the impedance was 
not accurately known. This study deals with impedance meas- 
urements using a Helmholtz resonator as a calculable stan- 
dard to determine the absolute accuracy and relative 
precision of the two techniques. A Helmholtz resonator was 
chosen because its various acoustical parameters could be 
independently measured or calculated and then compared to 
the experimental results. 

A. THEORETICAL BACKGROUND 

The term "acoustic impedance" can have different mean- 
ings for different situations. It is necessary to define 
exactly the different impedances involved in this study. 
This is necessary to prevent confusion later when the rela- 
tionship between the different impedances is used to deter- 
mine the calculated acoustical properties of a Helmholtz 
resonator. Normally this is not necessary since the specific 
acoustic impedance and the acoustic impedance differ by a 
constant wirh the dimensions of area. 

• Specif ic Aco usti c I moe dance 

The specific accustic impedance is defined as the 
ratio of acoustic pressure in a medium to the associated 
Eulerian particle speed. [Ref. 17] 

Z 5 = p/u ( 2 . 1 ) 

The MKS unit of specific acoustic impedance is ?a-s/m. 
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The specific acoustic impedance of a medium is a 
raai quantify for plane waves. However, this net rrue for 
the standing waves generated in an impedance rube. For ohis 
case, nhe specific acousoic impedance will be complex ana 
can expressed in the form 

Zi = r + jx (2.2) 

where "r" is the specific acoustic resistance and ”x" is rhe 
specific acoustic reactance for the medium concerned. 

2- ncoustic Impedan ce 

The acoustic impedance of a medium acting cn a 
surface cf area "A" is the ratio of the acoustic pressure at 
the surface to the volume velocity at the surface [Ref. 17] 

Zo. = P/U (2.3) 

where U = uA if u is normal to the surface of area "A". The 

HKS unit of acoustic impedance is Pa-s/ra^, which is also 
termed an acoustic chm. The acousric impedance and the 
specific acoustic impedance are related at the surface of 
area "A" by 

Za. = Z^/A (2.4) 

3 . Hec nan leal I mpedan ce 

The mechanical impedance is the ratio of the driving 
force to the resultant particle speed at the point where the 
force is applied [Ref. 17] 

z^ = F/u (2.5) 

The MKS unit of mechanical impedance is N-s/m, also referred 

to as a mechanical ohm. 
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4 . Hel mholtz Re s en ator 

A Helmholtz resonator can be treated as a iumpad 
parama'csr syscea. The rssona'cor consists of a vcluos "V 

with rigid walls and a neck of length "L” and cross- 
sectional area "A". The acoustic pressure in the volume 
provides the lumped stiffness ”s" if rhe wavelength is much 
greater than the cube root of the volume. The fluid in rhc'- 
neck provides the lumped mass as long as the wavelength is 
much greater than the length of the neck. Finally, the 

resistance is provided by the opening radiating sound, if 
the wavelength is much greater than rhe square root of the 
cross-sectional area of the neck. An additional resistance 
is introduced by viscous losses in the neck and irreversible 
thermal conduction to the volume boundries if the compres- 
sion of the gas is adiabatic and the walls remain isorh- 
ermal. However, since these resisrance sources are not 
easily calculated from "First Principles", they are 
neglected. The radiation resistance "R " can be calculated, 
however, its value is dependent upon the impedance of the 
medium which changes. To get the viscous loss resistance 
term the flow field in the neck is required. Therefore these 
resistance terms are cf little use in this calculation. 

The fluid in the neck has an affective mass given by 
[Ref. 17] 

m =/oAL' (2.6) 

where "L' " is the effective length of the neck and is a 

function of actual length "L" and the neck radius "a" and 
the flew field in the vicinity of the oriface. The effective 
mass is greater than the actual mass in the neck due the 
radiation-mass loading at the ends generated by this 
extended flow field. The stiffness of the resonator is 
given by [Ref. 17] 

s = (/ 3 C^) (A^ )/7 (2.7) 
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By driving rhe Helmholtz resona-or with a sound 
cf pressure amplirude ”P", the following equation 
inward displaos msnr "z” of the fluid in the nsck results 
[ 8ef . 17] 

sx = APe*'^*^ (2.8) 

This differsntal equation is the same as that for a driven 
oscillator and the irechanical impedance of the Helmholtz 
resonator can be expressed as [Ref. 17] 

z>„ = Hr + j (->m - s/o ) (2.9) 

where the mechanical reactance "X;„" is given by 

= com - s/to ( 2 . 10 ) 

Since the resonance occurs when X^ =0, the resonant 
frequency "uo" follows 

u>} = s/m ( 2 . 11 ) 

Using equations (2.6) and (2.7) 



vp/)/' 



The final result obtained from cancelling terms is 



- 




( 2 . 12 ) 



It should be noted that nowhere in the derivation of 
equation (2.12) is any restriction placed on the physical 
arrangement of the Helmholtz resonator. Figure 2. 1 shows the 
arrangement used for this study. As can be seen from 
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equation (2.12) ar.d the figure, the resonant frequency car- 
ts changed by changing nhs area of "neck" (diameter D, ), r.i-:- 

vciume of the resonator (change depth h) or the length of 

Recalling equation (2. 10) 






= u)m - s/f^ 



( 2 . 10 ) 

(2.11), th = 



and by factoring out "m" and using equation 
following result is obtained 

Since to = 2»7‘f the final form is obtained 
= 2rm(f - C/f) 

where f is the frequency in hertz. For convenience we define 
= f - f|/f (2. 13) 

so that 



X« = 27-mX:. 



(2. 14) 



It is now possible to derive the relationship 
between the specific acoustic reactance "Xj " and "X^,". Given 
a resonator as shown in Figure 2.2, divided into sections 1 
and 2 as indicated and recalling equation (2. 1) 



Zs = Pi /’i/ (2. 15) 

where p^ and u, represent the pressure and particle speed, 
respectively, in section 1 . This is the specific acoustic 
impedance of the medium in section 1. 

Next using equation (2.5) 






(2. 16) 



Converting Fj^ to a pressure at the boundary results in 
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(2. 17) 



= 



where p, and u^_ represent the pressure and particle sosed, 
respectively, in section 2. is the cross-sect ic r.al area 

of the resonator Oxeck. Since there is no physical boundary 
between sections 1 and 2, the pressures on both sides must 
be equal. 

Pi = Pz = P 

Also, the volume velocities must be equal across the 
boundry. Therefore 

u, = A^u^ 

u* = -^/-u, (2.18) 

Substituting equation (2.13) into equation (2.17) 
and regrouping the appropriate terms results in 







(2. 19) 



Finally, it can be seen that the quantity (p/u, ) is just the 
specific acoustic impedance and the result is 




Since the impedances are composed of the resistive and reac- 
tive parts, the reactive portions vary in the same manner. 
Therefore, 



= A.. 
K 



( 2 . 21 ) 
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It can be seen from eguatica (2.14) i:han 




2;7m- 






X 






It should be noted that the two techniques in this 
study provide the values for specific acoustic iapedar.ee 
normalized to the characteristic impedance "/jc". Therefore, 
the preceding equation becomes 



* w r 



and, finally, substituting for "m” from equation (2.6) and 
for "X_j/X,„ " from equation (2.21) 



^ ZIU'^ 



( 2 . 22 ) 



As can be seen, this equation does not depend on 

frequency and therefore remains constant as frequency 

changes. For convenience this constant is called "K^” with 

the subscript signifying the theoretical value for the ratio 
/ 

of "Xj//<?c]^" for a given Helmholtz resonator. 



K 



i-h ~ 






(2. 23) 



'• has been reduced to an expression which 
contains only directly measurable terms with the exception 
of "L"', the effective length of the neck. However, by using 
equation (2.12) "L*” can be determined from parameters which 

are measurable by separate means. Later in this chapter, 
"L’ '• will be looked at more closely for comparison with the 
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calculated "L"’ 
conf igurations , 



for 



ether 



Helmholtz 



resonator 



B. STANDING WAVS TUBS TECHNIQUE (SKI) 

The standing wave tube technique provides a convenient 
laboratory method for determining the absorption coefficient 
and specific acoustic impedance of materials. The method, 
however, is somewhat tedious when a large number of samples 
is measured or measurements of a single sample are required 
at a great number of frequencies. 



1 . Theor y 



The sample to be measured is placed at one end of 
the standing wave tube. At the other is mounted a loud- 
speaker. A microphone probe is inserted through the loud- 
speaker end and is attached to a microphone car. The car 
moves such that the microphone probe is moved along the axis 
of the tube. 

A pure tone signal is reflected from the sample. The 
reflected wave interferes with the incident wave and a 
standing wave is established. The microphone probe and car 
arrangement allows for the location and acoustic pressure of 
the nodes and anti-nodes to be measured. From the measure- 
ment of these amplitude and phase characteristics, calcula- 
tion of absorption coefficient and specific acoustic 
impedance of a sample is possible. 

From reference 6, the following are the equations 
needed for the determination of the desired acoustic 
parameters . 



X 



.TAn.‘£^l.LzTA~A ISil 

Ta,nh^ 



( 2 . 24 ) 



where 
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(2.25) 



'■P, = coth 



-I 



P, 



] 



and 





(2.26) 



and are the naximua and ainiriUn pressu 

standing wave tube and "d” is the distance to 
pressure nininum (node) from the sample. (NOTE: 
phone output voltage is the actual measurement 
only the ratio of the pressures is required) 



the first 
The micro- 
made since 



2. Eq u ipm ent Se t-u p 

Figure 2. 3 is a simplified representation of the 

equipment used in the standing wave tube technique. It 

includes the Bruel and Kjaer Standing Wave Apparatus Type 
4002, a HP 302A Wavs Analyzer and an amplifier and frequency 
counter. 

a. Bruel and Kjaer Type 4002 Standing Wave 

Appa ratus 

The apparatus is designed for easy and quick 
evaluation of acoustic materials by the standing wave 
method. A major advantage is that relatively small circular 
samples (10 cm. in diameter) are all that are required. 

The tube used had an interior diameter of 10 
cm., which covered a frequency range of 90 to 1800 hz. The 

lower limit was sat by the length of the tube and the upper 

limit determined by the cut-off frequency for the lowest 
non-plane wavemode of the tube. [Ref. 18] Two sample 

holders were also used with fixed depths of 1 in. and 2 in., 
respectively. The bases of the holders were very thick to 
present as rigid a termination as possible. 
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The micrcphone in the microphone car is a small 
dynamic microphone of adequate frequency range for the 
standing wave tube method. It is isolated by elasoic 
supports in the car ro reduce *: he efiecos of e.-ctsrnai 
airborne noise and impact sound or vibrarion. 

b. Hewlett-Packard HP 302A Wave Analyzer 

The Hewlett-Packard model 302A wave analyzer is 
a tunable voltmeter of high sensitivity over the frequency 
range of 20 Hz. to 50 kHz. The wave analyzer in the 3FO mode 
of operaticn becomes an audio signal generator combined with 
a tuned voltmeter. Further, since the tuned voltmeter has a 
very narrow pass band (6 Hz) , noise and harmonics of the 
signal have very little effect on the voltmeter reading. 

3 • Experi ments 1 Proce d ure s 

The standing wave tube technique was used in several 
different experiments. It was used with a rigid end as the 
sample to measure the "acoustic" end of rhe microphone probe 
tube and the sound speed. It was also used to determine the 
resonant frequency and reactive impedance of three Helmholtz 
resonators. Table I shows the physical dimensions of each of 
the resonators. 

a. Data Collection 

In all of the experiments the data was collected 
in the same manner. However, in the first experiment (rigid 
end) the only data collected was the location of the first 
two nodes and the frequency of the signal. Table II shows 
this data. 

For the experiments involving the Helmholtz 
resonators the values of the maximum and minimum microphone 
output voltages were obtained in addition to the frequency 
and node locations. This data was placed in Tables III-V 
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TABLE I 

H^lsholtz Resonator Physical Dimensions 



Resonator 


Volra ue 


SiSiS Diameter 


Neck Length 


no. 


cm . 


cm . 


cm. 


1 


195. 52 


2.52 


1.27 


2 


381 . 01 


2.52 


1.27 


3 


195. 52 


3.46 


1.00 



I 

I 

J 

b. Data Reduction 

In the first experiment, the measured distance 
between the two nodes was used to calculate the wavelength. 
From this, the quarter- wavelength was determined and 

compared to the distance from the sample to the first node 
as measured by the microphone car. Since the "sample” was a 
rigid end located 5.00 cm. from the measured end of the 
tube, the two distances should have differed by 5 cm. 
Additionally, the wavelength and the frequency were used to 
calculate the sound speed in the tube for each frequency 
measured. These results are- shown in Table II The average 
value of " d" suggests that the "acoustic center" of the 
probe tube is 1.7 + .5 mm. ahead of the physical end of the 
probe tube over the frequency range examined. 

For the experiments involving the Helmholtz 
resonators, the reduction of the data was somewhat 

different. The wavelength was determined by one of two 
methods depending upon whether two node locations were 
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TABLE II 








Standing «a 


vs Tubs Sound Speed 


T> A. 

L>a wci 




Freo . 


First Min. 


sec. Min. (;l^4) 




c 


Hz. 


c m. 


0 m • w til • 


cm . 


cm. /s 


299 


23.90 


81.74 28.93 


5.03 


34588 


400 


16.80 


59.95 21.58 


4.78 


34520 


600 


9.58 


38.38 14.40 


4. 82 


34560 


800 


5.96 


27. 53 10.78 


4.79 


34512 


1001 


3 .84 


21.03 8.60 


4.76 


344 14 


1 197 


2.19 


16.91 7.35 


5. 16 


35240 


1402 


1 .40 


13.64 6.13 


4.73 


34321 


1599 


0.54 


11.36 5.40 


4. 86 


34602 






average 


4. 87 


34595 






dsviat ion 


.150 


277 






per cent 


03.08 


0. 80 



measurable. In those cases where two or more node locations 
were measurable, the distance between them was used to 
calculate the wavelength. At the lower frequencies where 
only one node was measurable, the wavelength was determined 
from the frequency and the average sound speed found during 
the experiment. It should be noted that an attempt was made 
to use the location of the measurable anti-node to determine 
wavelength. However, this measurement was too inaccurate due 
to the nature of the waveform being measured. At the 
maximum, the values of amplitude are changing very slowly 
(second order) so that the measurement of the exact location 
of the maximum was very difficult. 3n the other hand, at the 
nodes, the values were changing more rapidly (first order) 
and the location could be easily and precisely determined. 
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TABLE III 



Standing Wava Tuba Heasurad Data - Resonator #1 



Freq. 


First 


Min . 


Sac . 


Min. 




Min. 


Maximum 


Hz. 


cm . 


uv . 


cm . 


U V. 


cm . 


U V. 


mv . 


250 


31.75 


65 










2.20 


350 


20. 85 


50 


69. 70 


60 






1.3 0 


400 


16. 65 


260 


59. 69 


260 






4.80 


450 


12. 3 0 


140 


50. 40 


130 


39.20 


175 


1 . 80 


500 


6. 70 


450 


41 . 20 


460 


76.30 


520 


2.30 


550 


1.70 


240 


30. 35 


290 


61 . 65 


190 


1 .60 


600 


. 55 


200 


22. 50 


180 


51 . 30 


135 


1.78 


650 


. 10 


200 


18. 42 


140 


45. 10 


165 


2. 60 


700 






15. 90 


34 


40.64 


31 


1.00 


750 






14. 20 


90 


37. 15 


1 00 


5.20 


1000 






9. 75 


30 


27. 15 


20 


3.40 


1 250 






7. 60 


5 


21 . 40 


5 


1 . 95 


1 500 






6. 75 


25 


17. 80 


5 


2. 50 






1 

I 



Standing Wa7e Tube Measured Data - Resonator #2 



Freg. 


First 


ain. 


Sec . 


Min. 


Third 


Min* 


Maximum 


Hz. 


cm . 


uv . 


cm . 


uv. 


cm . 


U V. 


mv. 


250 


27.88 


1750 










22. 1 


300 


19.69 


335 


76. 72 


330 






2.90 


350 


10.12 


1610 


59. 32 


1870 






8. 45 


375 


3.65 


1 280 


49. 63 


15 10 






6.70 


400 


1.59 


800 


40. 94 


915 


84 .65 


355 


4.35 


450 


.43 


480 


29. 90 


615 


68.58 


565 


8. 00 


500 


0.00 


1 64 


23. 96 


75 


58.64 


1 28 


2. 88 


750 


13.53 


24 


36. 35 


52 


59.50 


38 


2.9 1 


1 000 


9.75 


13 


26. 96 


15 


44.33 


39 


2. 43 


1 250 


7.44 


24 


21. 71 


43 


35. 12 


88 


2.00 


1500 


6.32 


30 


17. 93 


10 


29.37 


42 


4.42 



I I 
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TABLE V 



St 




Have Tub 


3 Heasured Da 


ta - Resosato 


r t 3 


Fr eg. 


Firs' 


t Hin. 


Sec . 


Hin. 


Third 


Min. 


Maximum 


Hz. 


cm . 


mv . 


cm • 


m V. 


cm . 


ra V. 


mv. 


200 


4 1.6 1 


3.2 0 










51.5 


250 


31.74 


29.5 










5 97 


300 


26.20 


2.90 


83. 95 


3. 15 






55.5 


350 


2 1.5 1 


3.90 


70. 04 


3.30 






65.5 


400 


18.42 


10. 0 


61 . 49 


9.40 






178 


450 


16.17 


1.52 


54. 44 


1. 70 






26 


500 


14.26 


4.05 


48. 19 


3.40 


82.21 


4. 05 


61 . 5 


550 


11.10 


5.40 


42. 92 


5. 30 


73.68 


5. 60 


68.5 


600 


9.18 


2.05 


37. 66 


1.98 


66.69 


1.70 


24.0 


650 


6.4 3 


6.05 


32. 79 


6.15 


60.39 


6. 10 


46.2 


700 


4.28 


4.00 


28. 62 


3.95 


52. 94 


3.80 


26.3 


750 


2.29 


2.90 


24. 54 


2.55 


47.49 


2. 64 


15.4 


800 


19.07 


3.95 


40. 41 


4.25 


62.37 


4. 15 


27.3 


850 


16.52 


1.10 


37. 32 


1.12 


57.42 


1. 10 


10.3 


900 


14.4 1 


3.10 


33. 30 


2.90 


52 . 80 


3.20 


47.5 


10 00 


1 1.67 


.80 


29. 00 


.90 


46. 18 


. 76 


22.5 


1250 


8.22 


.60 


21 . 86 


.45 


36.03 


. 52 


24.6 


1500 


6.53 


.65 


17. 96 


1.40 


29.49 


1 . 13 


60.5 



After the wavelength was determined, it was used 
to scale the distance to the first node (Table VI) . These 
values were plotted versus the frequency (Figures 2. 4-2. 6). 
The first two points above resonance were also plotted as 

d/A - .5 (2.27) 

A linear regression was then used with the two data points 
both above and below resonance to determine the resonant 
frequency ”f ". The resonant frequency for each resonator is 
also shown on Figures 2. 4-2. 6, 

Once the resonant frequency was obtained, 
was calculated for each frequency using equation (2.13). 

//3 c" was then calculated using equations (2.24-2.26). 
Finally, the ratio of "Xy//Oc" to "X^” (K«^p) was calculated 

for comparison to the theoretical value (K^^) determined from 
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TiBLE 71 



Standing Have Tube Data 



Free. 


resonator #1 


3 esonator #2 


5esoni 


Hz. 


d/j» 


<3 /A 


d, 


200 






.2416 


250 


.2437 


. 201 5 


. 2303 


300 


.220 1 


. 1726 


. 226 8 


350 


.21 34 


. 1028 


.2214 


375 




. 0397 




400 


.1934 


. 4683 (-.03 17) 


.2138 


450 


. 1585 


. 3865(-. 11 35) 


.2113 


500 


.0954 


. 3454 


.2099 


550 


.4848 (-.0 152) 




. 1774 


600 


.3906 (-. 1 094 




. 1596 


650 


.3452 




.1192 


700 


.3213 




. 0880 


750 


.3094 


. 2943 


.0507 


800 






. 4404 


850 






. 4039 


900 






.3754 


1000 


.2802 


. 2820 


. 3382 


1250 


.2754 


. 2688 


. 2950 


1 500 


.2706 


. 2751 


. 2844 



(-. 0596) 
(-. 096 1) 



equation (2.23). Tables VII-IX list the experimental values 
for each Helmholtz resonator. 

c. Error Analysis 

Before an analysis of the errors could be 
attempted, it was necessary to identify all of the sources 
of the errors. The errors were a result of the measurements 
of the voltage ar.d P^,^) , frequency and distances to the 

nodes . 

Hecalling equation (2.25) 

V' = coth"' (-6^-^--) (2.25) 
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TABLE 7II 

Standing Wave Tube Reduced Data - Resonator #1 







542.7 


Hz . 




Frequency 


c 




U 




Hz . 


cm. /s 


Pa-s/m 


S ^ X 


1 0-3 


250 


32050« 


-7.39 


-428. 1 


7. 


350 


34 1 95 


-4. 15 


-491.5 


8. 44 


400 


3443 2 


-2. 63 


-336.3 


7.82 


450 


34 60 5 


-1.52 


-204.5 


7.43 


500 


34800 


-.658 


-39.05 


7. 39 


550 


34430 


.094 


14.50 


6.48* 


600 


34 56 0 


.807 


109. 1 


7.40 


650 


34 6 84 


1. 45 


196.9 


7. 36 


700 


3463 6 


2. 07 


279.3 


7.41 


750 


3442 5 


2.55 


357.3 


7. 14 


1000 


34 80 0 


5. 21 


705.5 


7.38 


1250 


3450 0 


6.22 


1014 


6.13 


1500 


34650 


7. 70 


1304 


5.90 


average 


34560 




average 


7.57 


deviation 


174 




deviation 


. 39 


per cent 


.50 




per cent 


5. 15 


* Values not 


used for 


the calculation of th 


e averag 



it can be seen that for the values of ” 

involved, SWR is small and therefore, an error in the 

p 

standing wave ratio (SWR = ) has little effecr on rhe 

overall calculations. 

Looking at equation (2.24) 



- /2 24 ) 

and taking into account that ” " is small, equation (2.24) 

then reduces to 

^ / 
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TABLE VIII 



Standing 


Wave Tube 


Seduced Data - Resonator #2 






387.8 


Hz . 




Freguer.cy 


c 








Hz. 


cm . /s. 


Pa-s/m 


S~1 X 


10-3 


200 


34392* 


-3.86 


-551,9 


6.99 


250 


34 4 3 0* 


-2.43 


-351 . 6 


6.91 


300 


34542 


-1.46 


-201.3 


7.25 


350 


3444 7 


-.563 


-79.69 


7.06 


375 


33563* 


-.070 


-26.04 


2.69* 


400 


34 29 6 


.2 58 


24.03 


10. 74* 


450 


3351 6 


.745 


1 15. 8 


6.43 


500 


34620 


1.55 


199.2 


7.78 


550 


34557 


2. 08 


276.6 


7.52 


600 


3439 8 


2. 35 


34 9.4 


6.73 


650 


34437 


2. 80 


418.6 


6.69 


700 


3450 3 


3. 23 


485.2 


6.66 


750 


3450 8 


3.72 


54 9.5 


6.77 


1000 


34 56 0 


5.49 


84 9.6 


6.46 


1250 


3491 3 


8.51 


1130 


7. 53* 


1500 


344 70 


7. 1 1 


1400 


5.08* 


average 


34 44 4 




avera ge 


6.94 


deviation 


314 




deviation 


.41 


per cent 


. 9 1 




per cent 


5,9 


♦Values not 


used for t 


he calculat 


ion of the 


aver age 



( j 



or 









cot 'f’l. 



away fiom resonance. 

On the basis of this result, it was determined 
that the best data was obtained in the regions just away 
from resonance for the following reasons. First, since tan- 
goes to zero at resonance, ” terms begin to dominate and 
the errors in "SK” become large when compared to 
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TABLE IX 



Standing Wave Tube Reduced Data - Resonator #3 



Frequency 


c 


f^= 773.0 


Hz . 




Hz. 


cm. /s. 


Pa-s/m 


S“^ X 


10-3 


200 




-7.86 


-2788 


2.81 


250 




-6.93 


-2140 


3.24 


300 


34 65 0 


-5. 72 


- 1692 


3.38 


350 


33 99 9 


-5. 03 


-1357 


3.71 


400 


34 45 6 


-4.09 


-1094 


3.74 


450 


34 44 3 


-3. 78 


- 878 


4.31 


500 


33 97 5 


-3.63 


- 695 


5. 22* 


550 


34 41 9 


-2.07 


- 536 


3. 86 


600 


34506 


-1 . 50 


- 396 


3.79 


650 


35074 


-.8 70 


- 269 


3.23 


700 


34062 


-.600 


- 154 


3.90 


750 


33900 


-.3 13 


- 46.7 


6.70* 


800 


3464 0 


.427 


53. 1 


8.04* 


850 


34 76 5 


.6 16 


147 


4. 19 


900 


34 55 1 


1.06 


236 


4. 49 


1000 


3451 0 


1.58 


402 


3.93 


1250 


34 82 5 


3. 97 


772 


5. 14 


1500 


3444 0 


4.69 


1102 


4.26 


average 


3445 1 




average 


3.74 


deviation 


327 




deviation 


. 47 


per cent 


.95 




per cent 


12.6 


^Values not 


used for 


the calcula 


tion of the 


average 






Secondly, as the frequency gets further from resonance the 
cotangent function gets larger and the errors also get 
larger. It was for these reasons that some of the data in 
Tables VII-IX were deleted for the calculation of the 
aver ages . 



This completes 
wave tube technique. The 
later in this chapter along 
dual-microphone technique, 
technique needs addressing. 



the discussion of the standing 
final results will be examined 
with those for the Chung-Blaser 
One additional experimental 
The resonant frequency of each 
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Helmhcltz resonator was measured independently of 
method. This was done using a signal generator and 



to excixs the 


r eson 


atcrs. Using a 


small 


m i cr 0 p h o n 


resonators ~o 


detect 


xhe or assure cr 


1 a n g e 3 , 


the siqz-\ 


at or was tun 


ed to 


the frequency 


which 


resulrei 



greatesr pressure in the resonators. These 
frequencies were 549.5 Hz. for resonator 1, 398.5 

resonator 2 and 795.0 Hz, for resonator 3. 




in xhe 
resonant 
Hz. for 



C. CHONG-BLASER DUAL-MICROPHONE TECHNIQUE 

The dual-microphone technique developed by Chung and 
Blaser provides a quick method for measuring the acousxic 
impedance and reflection coefficient in rhe laborarory. 
However, as will be shown later, the technique suffers some 
degredaxion in accuracy compared with the results obtained 
by the SWT method. 

1 . The ory 

Chung and Blaser have shown [Ref. 1] that the 
complex reflection ccefficients at the two microphone loca- 
tions can be expressed in terms of the transfer functions of 
the incident and reflecxed waves. This expression is 



Ml . 



(2.28) 



where R^ (f) and R^ (f) are xhe complex reflection coeffi- 
cients evaluated at the first and second microphone loca- 
tions, respectively. H^* ( f ) and H,.(f) are the acoustic 
transfer functions associated with the incident and 
reflected wave components, respectively, evaluated at the 
two microphone locations. 
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For plane wave propagation (no mean flew) and 
neglecting wall losses, ( f ) and can b“= expr = ss-^d as 

(f) = (2.29) 

and 

(f) = (2. 30) 

where "k" is the wavenumber and "s" is the microphone 
spacing. 

Application of equation (2.28) resulrs in the 
reflection coefficient '*R'* on the surface of a sample 
material not at a microphone location being expressed as 





(2.31) 



wher e 
H 




(2. 32) 



and 

h; = e-*^' (2.33) 

and "1" is the distance from the firsr microphone to the 
surface of the sample material. Using equations (2.31-2.33), 
the reflection coefficient at the sample is 



R 



R, 



(2.34) 



where [Bef. 1] 



7/rG) - 



(2.35) 



H,^(f) is the acoustic transfer function between microphones 
1 and 2 and (f) and H,.(f) are the transfer functions of 
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the incident end reflected waves, respectively, between the 
two uiicrcphcne locations. 

fjsir.g the standard relation between reflection coef- 
ficient and ccnplex specific acoustic iffipedancs 



.^ 1 . = L.LtJ^X 

Ci~R) 



(2.36) 



and the relationships in equations (2.34) and (2.35), the 
complex specific acoustic impedance can be determined at a 
distance "1" from the first microphone [Ref. 1] 



-is.. = 



(2.37) 



C6s^£H /-5)3 A^/-s)3] t lt^J'’CosKi^i) 



H^^(f) can be expressed as 



(2.39) 

where C,^ and Q,, are the real and imaginary parts of the 
cross-power spectrum of the signals at microphones 1 and 2. 
S ,, (f) is the auto-power spectrum of the signals at 
microphone 1. 

Finally, Seybert and Ross show [Ref. 15] that the 
phase change upon reflection is given by 



0{t) 



tan~i 



C/rCO 



(2.4 0) 
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where Q^> ard the imaginary and real parts of the 

cross-power spectrum of the incident to reflected waves. 

2 • Equ iom en t Se t-u p 

Figure 2.7 shows a simplified diagram of the 
Chung-Blaser dual-microphone technigue equipment set-up. At 
the canter of the system are the Hewlett-Packard HF-35 
computer and HP-5420 Signal Analyzer System. 

a. The Hewlett-Packard HP-85 Computer 

The HP-85 computer is an 3 bit microprocessor 
which uses the BASIC computer language. It has 16K bytes of 
read/write memory of which 14,579 bytes are available for 
use by the operator. The unit used in the study had its 
memory expanded to 32K bytes by a memory module. The unit 
has a mass storage capability in the form of a magnetic tape 
reader/rscorder- It also includes a 127 millimeter diagonal 
black and white electromagnetic deflection CRT and a 32 
character per line thermal printer/plotter. To interface 
with other equipment, an I/O ROM and interface card were 
added to provide HP-IE (IEEE-488) interface capability. 

b. The Hewlett-Packard HP-5420 Signal Analyzer 

The Hewlett-Packard HP-5420A Signal Analyzer 
provides digital time-domain and frequency-domain analysis 
of low frequency analog signals. The HP-5420A system 
includes the HP-5447B Digital Filter Module, HP-544 10A 
Analog no Digital Converter, HP-5443A Keyboard and Control 
Module and the 544 1 A Display Module. The equipment used in 
this study also included an HP-10920A HP-IB interface acces- 
sory. This allowed its use with an HP-IB controller. It 
allowed the transfer of data and measurement set-up informa- 
tion to pass to and from the controller and allowed the 
reception of remote control commands. Using the controller. 
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the data was modified, processed, saved and recalled 
resal'^'ing in a high degree of flexibility. 

The HP-5U20A Signal Analyzer psrformes a varrety 
of ri me- do main and frequency-domain measurements. Ail meas- 
urements are done using digital signal processing. This 
involves filtering and sampling the input waveforms. The 
analog input is first filtered by a lew-pass filter and then 
passed re an analog to digital converter. Here, the waveform 
is sampled which results in the conversion of a voltage (at a 
csrrain point in time) into a numerical representation. 
These numbers are then filtered and processed digitally to 
produce the various measurements preformed by the analyzer. 

c. The Impedance Measurement Tube 

The tube used in this part cf the study was a 
piece cf 3 inch inside diameter PVC pipe 12 inches long with 
a 3 inch loudspeaker mounted at one end. Two holes were 
drilled in rhe tube and fitted with 7/8 inch Tygon rubing to 
hold the two microphones. The first was located 11.42 cm. 
from the sample end (opposite the loudspeaker) and the 
second was located 4.59 cm. from the first. The distances 
were measured to the center of the microphones. 

The microphones used were General Radio 1962 1/2 
inch elect ret-condenser microphones. These microphones have 
a sensitivity level cf -40 dB re 1 v/Pa and have rypically 
flat response curves from 20 Hz. out to 1800 Hz. which cover 
the frequency range of the study. 

3 • Ex p eri m e nt a 1 Pr oce d ure s 

The dual-microphone technique was used to determine 
the resonant frequency and the reactive impedance of the 
same three Helmholtz resonators used with the standing wave 
tube technique. A major difference between the two techni- 
ques made itself known almost immediately. The time to 
collect and reduce the data was vastly different. 
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a. Hicrophons Calibration 

The two niicroph on-r/pream plif iar channels had nc 
be calibrated for phase and aaplinude differences crier to 
use for the actual measurements. The cechnique used for this 
calibrarion was to place both microphones in idencical sound 
fields and record any frequency-dependent phase and ampli- 
tude differences. This calibration data was then used by rhe 
computer to condition che daca prior to the calculacicn of 
reactive impedance, phase angle and complex reflection coef- 
ficient. In this manner, the effect of differences between 
the two microphone channels on the data was minimized. 
Additionally, this method also allowed for imbalances in the 
A/D converter channels of the signal processor. 

The identical sound fields were obtained by 
mounting the two microphones at the same radial position at 
the sample end of the tube. This was based on che assumption 
that the sound field is radially symmetric at frequencies 
below the sloshing mode. This frequency is given by 
[Ref. 18] 

cx,c/d (2.41) 

where " is the solution to 

jLlmXu<d ^ 0 
dc<. 

and ”d" is the diameter of the 
corresponds to which is equal 

a cut-off frequency of 2650 Hz. 
used is given in Appendix A. 

b. Validation of Tubs Performance 

After calibration of the microphone channels, 
the microphones were placed in the normal positions in the 
wall of the tube and the proper operation of the system was 



(2.42) 

tube. The sloshing mode 
to .5861. This results in 
The calibration program 
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verified. k 4.03 inch long 
diameter, provided with a r 
a sample impedance. The ex 
cosfficienz of 1 and a knew 



piece of ?vc 
igid reflective 
psetsd results 
n phase change 



pips with the sam- 
end, was used a~ 
were a reflection 
equal to 




(2.43) 



The results are shown in Figures 2.3 and 2.9 Figure 2.8 
shows that the reflection coefficient is approximately 1 for 
frequencies below 2.5 kHz. except at 1.5 kHz., were the 
phase goes through 0 degrees yielding an indeterminancy 
which is discussed in section e. which follows. Above 2.5 
kHz. non-plane wave modes degrade the results. The measured 
slope of the phase shift versus frequency shown in Figure X 
is -53 degrees ± 1 degree. Using equation (2.41), this slope 
produces a sound speed of 356 + 1.34 m/sec, which is consis- 
tant but significantly higher than the average of the four 
values obtained by the SWT technique listed in Tables 
II,VII-IX (c = 345. 1 2 ± .38 m/sec). 

c. Data Collection 

Each of the three Helmholtz resonators was 
placed at the sample end and the data collected. The data 
was tabulated in Tables X-XII. The collection and recording 
of the data was done completely by the computer and the 
program is given in Appendix B. 

d. Data Reduction 

After the raw data was collected the first step 
was to determine the resonant frequency of the Helmholtz 
resonator. The data collected allows this to be done in two 
different ways. Both were used. First, the specific reactive 
impedance "X^ //>c" was plotted versus the frequency and a 
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TABLE X 

Dual-Microphone Measured Data 



Resonator #1 Run 1 



’T 



! 



I 

I 

? 

I 



f 4 = 58 7. 3 ± 2.2 Hz, 



Phase An^le 
deg. 







Hz. 


Pa-s/m 



250 


-5. 87 


350 


-2. 83 


4 00 


-1 . 91 


450 


- 1. 23 


500 


-.671* 


550 


- . 2 56* 


600 


. 1 14* 


650 


. 395* 


700 


.6 76 


750 


.9 19 


1 000 


1 . 86 


1250 


2. 66 


1 500 


3. 03 



Run 2 



Fr eg . 




Hz. 


Pa -s/m 



250 


2. 13 


300 


-1 . 58 


350 


-2. 14 


400 


- 1 . 64 


450 


- 1. 07 


500 


-.591 


550 


-.2 24 


6 00 


. 076 


650 


. 3 50 


700 


.6 29 


750 


.8 90 


800 


1 . 19 


900 


1. 52 


1 000 


2. 22 


1 250 


2. 96 


1 500 


3. 17 



0 

- 23.9 

- 41.7 

- 65.6 
-105.2 
-147.6 

166. 2 
135.6 
1 10. 7 
93.9 
55, 2 
43. 2 
36. 5 



f, = 59 0. 6 ± .2 



Phase Angle 
deg. 



19. 2 

- 8. 4 

- 27.4 

- 52.2 

- 78.6 
-113. 5 
-152. 6 

170. 8 
140. 3 
114. 6 
96.0 

79.9 
66.5 
47. 8 
37. 1 

34.9 



36 0 + Pha se Ana. 

deg. 



360 
336. 1 

318. 3 

294. 4 
254.8 
212. 4 



Hz. 



360 + Phase An_g. 

deg. 



379. 2 
35 1 . 6 
332. 6 
307. 8 

281. 4 
246. 5 

207. 4 
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regression 


was used 


d 


ids cf the 


origin (4 




edarce was 


equal to 




y since ane 


reacr.v g 




resonance. 


The graphs 



W h 9 ~ <r 



;ir. ts on 

V - 



or 



each resonator are shewn in Figures 2.10-2.14 

The second method for determining the resonanr 
frequency involved the use of the phase angle between the 
incident and reflected waves. When *he phase angle was 
plotted versus frequency (Figures 2.15-2.19) resonance was 
obtained from the graph at the point where the phase angle 
was 180 degrees. Again a linear regression was used with the 
2 points on either side of the 180 degree phase angle. 
Figures 2.15-2.19 also show the frequency determined. The 
resonant frequency used in subsequent calculations was the 
average of the two frequencies determined from the two 
methods. The reproducibility of the results can be estimated 
from the ratio of the standard deviation to the average for 
the four (two) separate raeasuremenrs of the resonant 
frequencies for resonators 1 and 2 , which were 0.5® and 
2 . 0 %, respectively. 

Once the resonant frequency was determined, it 
was used along with equation (2.13) 



X * 



= f - f.Vf 



(2. 13) 



to obtain "XJi". Using this along with the reactive impedance 
ratio "Xj/pc", the experimental value for the rario 



K = 



(2. 42) 
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TABLE XI 



Dual-a 


icrophone 


Measured Data - 


Resonato 


r #2 Run 1 






435.8 i 


3.8 Hz. 








Pha^ Angle 


360 


+ Phase An^ 


Hz. 


Pa-s/n» 


deg. 




deg. 


200 


. 2 15 


7.5 




367. 5 


250 


-4. 24 


- 1.6 




358. 4 


300 


- . 1 51 


- 8. 3 




351 . 7 


350 


-.4 36 


- 47.3 




312. 7 


400 


- .206 


-121.7 




238. 3 


450 


. 1 00 


164. 6 






500 


.3 86 


133.0 






550 


.667 


109.7 






600 


. 893 


94.8 






650 


1 . 07 


85.0 






700 


1. 27 


76. 0 






750 


1 . 40 


70.4 






800 


1. 60 


63.8 






900 


1 . 84 


56.7 






1000 


2. 12 


50. 1 






1250 


2. 57 


42.2 






1500 


2. 83 


38. 1 






Run 2 




f . = 421.7 + 


. 3 Hz. 




Frj2* 




Phase Angle 


360 


+ Phase Ang 


Hz- 


?a-s/m 


deg. 




deg . 


200 


. 1 63 


1.6 




361. 6 


250 


-2 . 03 


- 41.7 




318. 3 


300 


-.9 87 


- 87.4 




272. 6 


350 


-.4 34 


-132. 5 




227. 5 


400 


-.082 


-170.5 




189. 5 


450 


. 1 71 


160.0 






500 


. 4 12 


134.4 






550 


.5 93 


116. 9 






600 


. 8 13 


100. 7 






650 


.9 93 


89.8 






700 


1. 27 


75.7 






7 50 


1 . 50 


67. 1 






800 


1 . 78 


58. 4 






9 00 


1 . 91 


54.8 






1000 


2. 65 


4 1.2 






1 2 50 


3.35 


33.2 






1 5 00 


3. 38 


33. 0 
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TABLE XII 

Dual-Iiicropbcne Haasured Data - Resonator 53 



f, = 94 7. 1 * 1.0 Hz. 



Fre^. 


in 


Pha^ Anole 


360 + Phase 


Hz. 


?a -s/m 


deg. 


deg. 


400 


. 7 19 


22. 1 


382. 1 


450 


. 1 74 


5.6 


365. 6 


500 


-. 039 


. 9 


359. 1 


550 


-.5 97 


- 15.4 


344. 6 


600 


-1. 01 


- 25.4 


334. 6 


650 


-1 . 06 


- 34.6 


325. 4 


700 


-.921 


- 56.4 


303. 6 


750 


- . 8 34 


- 69.4 


290. 6 


800 


-.6 25 


- 94.6 


265. 4 


850 


-.3 80 


-126.9 


233. 1 


900 


-. 201 


-150.7 


209. 3 


950 


.0 48 


173. 5 




1000 


. 1 96 


155.2 




1050 


.3 85 


135.5 




1 100 


.571 


118. 1 




1200 


1 . 05 


86. 5 




1300 


1 . 34 


73.0 




14 00 


1 . 57 


64.7 




1 500 


1. 69 


6 1.2 








was calculated for comparison with the theoretical value 
obtained from equation (2.23) 






(2.23) 



Tables XIII-XV show the results for for each of the 

resonators tested. 

e. Error Analysis 

The Chung-Blaser dual-microphone technique had 
some inherent errors built into it by virtue of the equip- 
ment construction and mathematical techniques. First, the 
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TABLE XIII 



Dual-aicrophone Reduced Data - Resonator #1 Ran 1 



Frequency 


U 






Hz. 


S~ 1 




X 10- 


250 


-11 42 




5 . 14« 


350 


- 644.2 




4. 39 


400 


- 470.0 




4.06 


450 


- 323.3 




3.80 


500 


- 196.0 




3.42 


550 


- 82.7 




3. 10* 


600 


20 . 0 




5.69* 


650 


114.6 




3.45 


700 


202.9 




3.33 


750 


286.0 




3.21 


1000 


6 52.0 




2.85 


1250 


971.6 




2.74 


1500 


12 68 




2.39 






av erage 


3. 52 






deviation 


. 79 






per cent 


22.4 


Run 2 








Frequency 


Xi 






Hz . 


s- 1 




X 10- 


250 


-1146 




1 . 86* 


300 


- 863.1 




1. 83* 


350 


- 646.9 




3.31 


400 


- 472.3 




3. 47 


450 


- 325.4 




3. 29 


500 


- 197.9 




2. 97 


550 


- 84.4 




2.65* 


600 


18.5 




4.12* 


650 


113.2 




3.09 


700 


201 . 5 




3. 12 


750 


284.8 




3.13 


800 


363.8 




3. 27 


900 


512.3 




2. 97 


1000 


651 . 1 




3. 41 


1200 


909.2 




3. 26 


1400 


1151 




2.75 






average 


3. 17 






deviation 


.21 






per cent 


6. 62 



♦Values not used for the calculation of the average. 
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TABLE XIV 

Dual-iJicrophone Beduced Data - Resonator #2 Run 1 



Frequency 








Hz . 


s~ ^ 




X 10-3 


200 


- 749.6 




. 29 * 


250 


- 509.7 




8 . 32 * 


300 


- 333.1 




, U 5 * 


350 


- 192.6 




2 . 26 


400 


- 74.8 




2.75 


450 


28.0 




3 . 58 * 


500 


120.2 




3.21 


550 


204 . 7 




3 . 26 


600 


283 . 5 




3.15 


65 0 


357 . 8 




2.99 


700 


428.7 




2.96 


750 


496 . 8 




2. 82 


800 


562 . 6 




2 . 84 


900 


689.0 




2.67 


1000 


810.1 




2.62 


1200 


1042 




2.47 


1400 


1264 




2 . 24 






av erage 


2.79 






de viarion 


. 33 






per cent 


11.8 


Run 2 








Frequency 








Hz. 


s- ‘ 




X 10-3 


200 


- 689.6 




. 24 * 


250 


- 46 1.7 




4 . 40 * 


300 


- 293.1 




3 . 37 


350 


- 158.3 




2 . 74 


400 


- 44.8 




1 . 83 * 


450 


54.6 




3 . 13 


500 


144 . 2 




2 . 86 


550 


226 . 5 




2 . 62 


600 


303 . 5 




2 . 69 


650 


376 . 3 




2 . 64 


700 


445 . 8 




2 . 85 


750 


512 . 8 




2 . 93 


800 


577.6 




3.08 


900 


702 . 3 




2 . 72 


1000 


822 . 1 




3 . 22 


1200 


1052 




3.18 


1400 


1273 




2 . 66 






average 


2 . 91 






deviation 


. 25 






per cent 


8.59 



♦Values not used in nhe calculation of the average. 
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TABLE XV i 

Dual-Hicrophcce Reduced Data - Resonator *3 I 



Fr equency 


Xi, 






Hz. 


s- 1 




X 10-3 


400 


- 1844 




. 39 * 


450 


- 1545 




» ♦ 


500 


- 1295 




. 03 * 


550 


- 1082 




. 06 * 


600 


- 895.9 




1 . 13 


650 


- 730.9 




1 .45 


700 


- 582.2 




1 . 58 


750 


- 446.8 




1 . 87 


800 


- 322.0 




1.94 


850 


- 206.0 




1 . 84 


900 


- 97.3 




2.07 


950 


5 . 2 




9 . 24 * 


1000 


102 . 4 




1.91 


1050 


195 . 2 




1 . 97 


1100 


284.0 




2.01 


1200 


452.0 




2 . 32 


1300 


609 . 6 




2.20 


14 0 0 


758.9 




2.07 


1500 


90 1.6 




1 . 87 






av erage 


1.91 






deviation 


. 23 






oer cent 


12.0 



* Values not used in the calculation of the average. 



construction of rhe tube limits the range of frequencies 
over which good dara could be obtained. The diameter of the 
tube determined the upper limit of measurable frequency. 
Since the tube was valid only for a plane wave, the frequen- 
cies measured had tc be below the curoff frequency for rhe 
tube. For the tube used in this study the cutoff frequency 
was approximately 2600 Hz. The construction of the rube also 
limits the effectiveness of measurements at the lower 
frequencies. The critical parameter in this case is the 
spacing between the microphones. For a given separation as 
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ths wavelength gets longer (frequency lower) , the differ- 
ences in the phase and ainplitudes at the two nicr''th.rn‘= 
locations becaae small and the noise dominated. The low 
frequency limit for the tube was about 200 Hz. 

Finally, the mathematical technique resulted in 
an indeterminacy in the reflection coefficient. This 
occurred when the phase angle was equal to zero. >Jhen this 
occurred equation (2.35) reduced to 

R = 0/0 

which is indet er oina nt. The result was that data obtained 

for the frequencies about 50 Hz on either side of the point 
where the phase angle equalled zero was greatly in error. 

D. CCMPAEISON OF THE TWO TECHNIQOES 

The two techniques were consistent with each ether and 
with the expected theoretical values. First, the techniques 
will be locked at fer comparison with the the theoretical 
results and then for cemparisen with each other. 

• Sta ndi ng Wave Tube Techn ique 

Figures 2.20-2.22 show the comparison of the experi- 
mental and theoretical values for ”Xj4<jc''. The theoretical 
values were obtained using (eqn. 2.23) and the calcu- 

lated values for "X^" (Tables VII-IX) . These were tabulated 
in Tables VII-IX As can be seen in Figures 2.20-2.22 the 
values compare favorably with a divergence noted at the 
higher frequencies away from resonance. This divergence was 
explained in the error analysis in section B.3.c. 

Table XVI shews the resonant frequency and both 
constants for each resonator, along with the resonant 
frequency determined by the independent method discussed in 
section B.3.c. As can be seen, were within 
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the calculated error for each. The average per cert error 
between the two values for ”K '* was found to be 4.13?; fcr 
the standing wave tube technique. The average per cent error 
for the resonant frequency was found to be -2.23%. 

2 • Chu nq- 3 las er Dual- W ic r ophon e Techn^gue 

Figures 2. 23-2.27 show the coniparison of rhe experi- 
mental and theoretical values for The th=^oretical 

values were obtained in the same manner as in the SWT tech- 
nique. The graphs shew a good agreement except at the lower 
and higher frequencies away from resonance for reasons 
previously noted in section C.3.e. 

Table XVI shows the resonance frequency and both 
values obtained for the dual-microphone technique. The 
values determined fcr and "K/.^ " do not in all cases 

overlap within the areas of uncertainty as was the case with 
the SWT technique. As a result, the average per cent error 
~ X 100 was found to be 18.4%. The average per 

cent error between the independennly measured resonann 
frequency and that measured by the dual-microphone method 
was found to be +11.4%. 

3. St and ing Wave Tu^ vs. Dua l- M ic r op h one 

The greatest discrepancy between the two menhods was 
the measured resonant frequencies. The resonant frequency 
measured by the dual-microphone nechnigue was always higher 
than the resonant frequencies measured by the other two 
methods. The average error of +11.4% was significantly 
larger than the error for the standing wave tube method of 
-2.23%. Additionally, the determination of tire constants 
for the dual-microphcne method produced larger discrepancies 
between the theoretical and experimental values. Finally, 
the calculated sound speed for the dual-microphone method 
was 3.2% higher than the SWT method which was in good 
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TABLE XVI 

Ccinparison of Results 



RES #J1 
549.5 i 1.5 



Frsq. 

Hz. 



Hz. 



;e < 

X 1 0-3 



Hi, 

sec. 



K 



ejtp 



sec. 
1 0-3 



fo 



( 

X 100 1 



SWT 


54 2. 7 


7.3 5 


+ . 18 


7.57 


± .39 


+ 2. 99 


Dual run 1 


589. 5 


3.6 7 


t .10 


3.52 


^ .79 


+ 4.09 


Dual run 2 


59 0. 5 


3.66 


+ .10 


3.17 


1 .21 


-13.4 


RES #2 














398. 5 ± 3.5 


Hz. 












SWT 


387. 4 


7.37 


± -IS 


6.94 


+ .41 


-5.83 


Dual run 1 


43 2. 0 


3.45 


± -33 


2.79 


± .33 


-19.1 


Dual run 2 


422. 0 


3.6 8 


.07 


2.91 


j, .25 


-20.9 


RES #3 














795.0 + 4.0 


Hz. 












SWT 


77 3. 0 


3.6 1 


t .09 


3.74 


+ .47 


+ 3.60 


Dual run 1 


97 4. 4 


1.4 2 


± .04 


1.91 


+ .23 


+ 34.5 



£JI CENT ER£0R in reson ant frequency 

R esen at or #1_ - 54 9.5 Hz . 

SWT - 542. 7 Hz . - 1.24% 

Dual - 590 . 3 Hz . ■►7.42% 

Re sonato r #2 ; 398.5 Hz . 

SWT - 387.4 Hz. -2.79% 

Dual - 428.8 Hz. *1,60% 

Res onat or #3 - 22£.0 Hz . 

SWT - 773. 0 Hz . - 2.77% 

Dual - 974.4 Hz. •►22.6% 
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agreement winh the expected sound speedy based 
ambient temperature. 



on the 
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WOODEN DISK 




Di D2 



Figure 2. 1 



Helmholtz Resonator Physical Set-Dp. 



Figure 2. 
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2 Theoretical Helmholtz Resonator Arrangement. 
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Figure 2.3 Standing Wave Tube Equipment Set-Op- 
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Figure 2.4 Resonator #1 Graph of d/^ vs. Frequency. 
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Figure 2.5 Besonator #2 Graph of vs. Frequency. 
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Figure 2.6 Resonator #3 Graph of d/^ vs. Frequenc 
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Figure 2.7 Dual--Microphone Technique Equipment Set-^up. 
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FREQUENCY-KHZ 

Figure 2.8 Reflection Coefficient vs. Frequency - Rigid End 
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Figure 2.9 Phase Angle vs. Freguency - Rigid End. 
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Figure 2.10 Dual-Microphone vs. Freq. - Resonator #1 Run 
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Figure 2.11 Dual-Microphone vs. Freg. - Resonator #1 run 
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Figure 2.13 Dual-Microphone "X^/oc'' vs. Freg. - Resonator #2 Run 
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Figure 2.14 Dual-Microphone "X^Xoc" vs. Freg. - Resonator #3. 
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Figure 2.15 Dual-Microphone Phase ftng. vs. Freg. - Resonator #1 Runi 
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FREQUENCY ~ KHZ 

Figure 2.16 Dual-Microphone Phase Ang. vs. Freg. - Resonator #1 Run 



320 + 



o 



Osl 



O 





4 1 1 1 h 



CD 


CD 


o o o 


O 


CD 


CXD 




CD LD CSJ 


CO 


cr 


OJ 


CNI 


CNI r— 1 rH 








03Q 


- 310NV aSVHd 







65 



Figure 2.17 Dual-Microphone Phase Ang. vs. Freg. - Resonator #2 Run 



360 + 



o 




^ 1 1 1 1 1 1 ^ 



CZD 


CD 


CZD 


o 


CD 


CO 


CD 


CD 


CNI 


CO 


cr 


cz: 


UD 


CNJ 


CO 


cr 


hO 


CNI 


CNI 


CNI 


rH 









93Q - HHONV 3SVHd 



66 



Figure 2.18 Dual-Microphone Phase Ang. vs. Freg. - Resonator #2 Run 
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Figure 2.19 Dual-Microphone Phase Ang. vs. Freg. - Resonator #3. 
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Figure 2.20 SWT Comparison of Reactive Impedances - Resonator #1. 
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Figure 2.21 SHT Cooparison of Reactive Impedances - Resonator #2. 
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Figure 2.22 SWT Comparison of Reactive Impedances - Resonator #3. 
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Figure 2.23 Dual-Mic. Comparison of Reactive Impedances 



< 




od/^X 



M 

X 



I 

> 

u 



UJ 

X 

o 

UJ 

cn 

u. 



CN 

a 

a 

a: 






12 



Figure 2.24 Dual-Mic. Comparison of Reactive Impedances - Res 
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Figure 2.25 Dual-Mic. Comparison of Reactive Impedances 
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Figure 2.26 Dual-Hic. Comparison of Reactive Impedances - Res. #2 Run 
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Figure 2.27 Dual— Mic. Comparison of Reactive Impedances - Res. 



III. CONCLOSIONS 



Throughout this study both techniques displayed advan- 
tages and disadvantages. The criteria for saying one is 
"better" than the other depends upon the purpose fcr which 
the information is needed. No attempt will be made to decide 
which technique is "better". The strengths and weaknesses 
will be pointed out and the potential users of this infcrma- 



tion can 


decide which technique more 


closely fits 


their 


needs. 








It was 


clear from the results in 


Table XVI that 


the 


Stan din g 


wave tube technique was the 


mors accurate 


when 



compared to the expected theoretical values for the 
Helmholtr "standards". In addition, the SWT technique had a 
higher average relative precision than the dual-micrcphcne 
technique. The average ‘fcr the SWT was 4.2^ while the 

average relative precision for the dual-microphone technique 
was 18.4%. We believe the reason for the factor of 

fcur-and-one-half better results for the SWT technique was 
due to one major characteristic of the two techniques. The 
standing wave established in the two tubes were theoreti- 
cally the same. However, the locations where the appropriate 
measurements were obtained were different for each tech- 
nique. Fcr the SWT technique the measurements were always 
taken at the optimum locations, ia. a maximum or a minimum 
pressure point. Therefore, it was possible to obtain the 
best data available. For the dual-microphone technique the 
measurements were taken at locations on the standing wavs 
predetermined by the microphone locations. These locations 
were only optimized when the microphone separation was equal 
to an odd multiple of a quarter-wavelength cf the frequen- 
cies measured. Therefore, the data obtained was not 
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nscsssarily the best available and the result was a degrada- 
tion in accuracy and relative precision for the iual- 
micropbone technique. This does not provide an ob'virus 
explanation for the error in measured sound speed or signrf- 
icant error in measured Helmholtz frequency obtained by the 
dual-microphone method. 

It would appear that the SWT method was the more desir- 
able based on accuracy and relative precision. This may not 
be the case. When numerous samples were involved or the 
results desired over a large frequency range, the SWT method 
became tedious and cumbersome. It involved the manual taking 
and reduction of data. The dual-microphone technique was a 
great deal faster and easier to use. The time to obtain the 
results for one sample over a frequency range of 200-1500 
Hz. was approximately 20 minutes for the dual-microphcns 
technique and 6 hours for the SWT technique. 

Thus when choosing which technique is more desirable for 
a given application, it must be decided whether speed or 
accuracy is more important. The degradation in the accuracy 
and relative precision for the dual-microphone technique 
might be considered minor in those applications where speed 
is important. 
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would allow the use of a na 
in a reduction in the noise 
useful redundancy as the wa 
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Although hare we only d 
would also be possible ro 
resistive component of the 
additional measurement. If 
nator is measured and the 
known, rhe absolute resist! 
can then be compared with t 
Finally, several chan 
dual-microphone technique 
First, the tube should be 
phone locations so that th 
could be varied somewhat t 
frequency range. Additiona 
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control, 
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he maximum 
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m noise, 
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rbow band 
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velength 
• 
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extract an 
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ve impedan 
he impedan 
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sary 
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0 

(D 
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e dual-micro 


use 


of a p 


It would be 
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use of 



study were 
to improve 
o adapt the 
be acccm- 
ive on the 
tically set 
values and 
tion. This 



phone tsch- 
ure tone as 
possible to 
a pure tone 
filter and could result 
nee and would include a 
and thus phase differed 

reactive impedance, it 
absolute value of the 
resonator with only one 
ty factor of the reso- 
reactive impedance is 
ce is fixed. This value 
ce measurement, 
be made to adapt the 
a water filled tube, 
ed with numerous hydro- 
on between hydrophones 
the data for a given 
sa of a bubble of known 
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APPENDIX A 

MICROPHONE CALIBRATION PROGRAM, 
option base 1 



I 


matrix 


1 r^l 


holds HP542 


:0 data 


bio 


ck header 






1 


ma L.rz-X 


’ A' 


holds auto 


an d cr 0 


33 


spect run 


;5 ii 

u. C. . 


^ a 


1 


matrix 


' C 


is call brat 


:ion vac 


tor 








t 


addrss 


s of 


HP-5420 


is 


704 










di 


.m W (2, 


16), 


A (3,256) , 


, C (2 


:,256) 










output 704 ; 


"RS" 3 S = 


= SPC 


ILL (704 


) 5) 


wait 500 


3 


s#o 


1060 


















f 


HP5420 


t es 


t setup: 














f 


60 St a 


ble 


averages 














f 


channe 


1 1 


active 














! 


tr ansf 


er f 


unction 














t 


ran don 


s ig 


na 1 














1 


1 . 6 kHz 


ban 


dw idth 














on 


itput 7 


04 ; 


"60,1 A7 


CH 


5 FR 2 


SG 


1.6 KV BW 


0 


CH 



10 RG" 

S = SPOLL (704) a wait 500 3 if S #) then 1140 
if s#0 then 1 l 40 
! start data collection 

disp "press end line to continue" a input A$ 
output 704 ;"ST" 

£=SPOLL (704) 3 wait 5 00 3 if S#68 then 1190 
disp "press end line to continue" 3 input A$ 
gosub 1330 3 gosub 1600 

! create data file'CALDAT* if not present on current 
tape 

on error goto 1240 3 create "CALDAT" , 1 , 4096 
off error 3 assign# 1 to "CALDAT" 
print# 1,1 ; C (,) 

! close fils to record data oi tape 
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1270 

1280 

1290 

1300 

1310 

1320 
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1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1580 

1590 

1600 

1610 



assigr.# 1 to 

disp "calibraticri data stored" 

disp "end program" 

clear 7 

stop 

and 

! collect data 

i transfer auto spectrum data to HF35 
output 704 ;'M TC 1 FM 1 AU" 

S = SP0LL (704) 2 wait 500 2 if S#4 then 1360 

output 704 ;"501SA" 

S=S?0LL (704) 2 wait 500 2 if S#32 then 1380 

for 1=1 to 16 

! enter 5420 header information into matrix W 
enter 704 ; W (1,1) 

next I 

! enter data in matrix A 

for 1 = 1 to 256 2 enter 704 ; A (1,1) 2 next i 

! transfer cross spectrum data to H?85 
output 704 ;"RS" 

S=SF0LL (704) 2 wait 500 2 if S#32 then 1470 

output 704 ; "CR" 

S=SP0LL (704) 2 wait 500 2 if S #0 then 1490 

output 704 ; "501SA" 

S = SPOLL (704) 2 wait 500 2 if S #3 2 then 1510 

for 1=1 to 16 2 enter 704 ; W(2,I) 2 next I 

for 1=1 to 256 

! enter real cross spectrum data to matrix ‘A’ 
enter 704 ; A (2,1) 

! enter imaginary cross spectrum data to matrix 'A' 
enter 704 ; A (3,1) 

next I 
return 

! calculate calibration vector 
print "calculating" 



8 1 



1620 for 1=1 to 256 
1630 C (1,1)=?. (2,1) /A (1,1) 
1640 C (2,1) =.A (3,1) /.A(1 ,1) 
165 0 next I 
1660 return 



% 

0 ? r«*z to 5 usar 

(7 • » » *S f f vSi A« (T »f > t> i 

t it9i^ oe^f 
• ii!^e3 o5if 



10 
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210 

220 

230 

240 

250 



AP^SNMX B 

DO AL-MICBOPHO >1E COMPOrSR PROGRAM 

laddrsss of HP5420 si 704, plottsr is 705 
! list of variables 



1 A(1,) 


mic 1 autospectrua 


• A{2,) 


cross spectrum, real part 


1 A(3,) 


cross spectrum, imaginary part 


t C(1,) 


calibration vecror, real parr 


! C(2,) 


calibration vector, imaginary part 


! RH 


real part of acoustical transfer function 


! IH 


imaginary part of acoustical transfer 


funct ion 




! A1() 


reflection coefficient 


! PHI 0 


phase change 


! R() 


impedance, real part 


! M() 


impedance, imaginary part 


option base 


1 3 dim A (3,256) ,A 1 (256) ,C (2,256) , M (256) , 



R (256) ,PHI (256) 

I all following air constants and calculations are in SI 
units 



13=295. 15 


! temperature 


G3=1. 4 


! specific heat ratio 


P3=101325 


! air pressure 


R3=237.0537 


! gas constant 



C=SQR (G3 *R 3*T3) I speed Of sound in air 
RO=P3/ (R3*T3) ! air density 

SOC=RO#C 

assign# 1 to "CALDAT'* a) read# 1,1 ; C(,) 



disp ’’enter 


number of materials to be measured” ; Sin put 


18 




for 19=1 to 


18 3 print ”we will now measure material #” 


;I9 
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26 0 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 



disp "if sample is ready, press end line" 3 input A3 
cutout 704 d S=SFOLL (704) 3 wait 500 3 if S#0 

’thsr. 27 0 

output 704 ;"30,1 AV 1 CH 5 FS 2 SG 1.6 KV 5W 0 CH 0 AC 
AC 10 RG" 

S=SPOLL (704) 3 wait 500 3 if S#0 then 290 
output 7 04 ;»ST MR" 

S=SFOLL (704) 3 wait 500 3 if S#68 then 310 
disp "end daxa collection" 
output 7 04 ;"1 FM 1 TC 1 AU" 

S=SPOLL (704) 3 wait 500 3 if S#4 then 340 
J=1 3 gosub 530 

output 704 ;"RS" 3 S=SPOLL (704) 3 wait 500 3 if S#0 

then 360 

output 704 ;"CR" 

S=SPOLL (704) 3 wait 500 3 if S#0 then 380 
gosub 530 

output 704 ;"RS" 3 S=SPOLL (704) 3 wait 500 3 if S^O 
then 400 

disp "off to calculate! this will take several 
minute s" 
gosub 600 

disp "press end line to store impedance data on tape" 
input A$ 

! the following statements save data on tape 
create "IMPDAT" ,256, 2 4 

assign# 1 to "IMPDAT" 3 for 1=1 to 256 3 F=(I-1)*D1 3 
print #1, I ; F, R (I) 3 next I 3 assign# 1 to "IMPDAT2" 
disp "end data collection for set";l9 
next 19 

disp "end program" 
and 

! start data collection 
output 704 ;"501SA" 

S=SPOLL (704) 3 wait 500 3 if S#32 then 540 
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550 

560 

570 

580 

590 

600 
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620 

630 

640 

650 

660 

670 

680 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 

810 

820 

830 



for 1 = 1 to 16 a) enter 704 ; 'li d) next I 
if J=2 then 580 

for 1 = 1 to 2 56 n) enrer 704 ; i (J,I1 ansxt I 3 orir.t 

"auto-spectral data transferred" 3 J=J+1 3 renurn 

for 1=1 to 2 56 a enter 7 04 ; A (2,1) 3 enter 704 ; A (3, 1) 

anext I a print "cross-spectral data transf erred" 3r eturn 

! calculate impedance 

disp "calculating" 

for 1=1 to 256 

C1=C (1 , 1)*C (1 ,1) +C (2,1) *C (2, 1) 

TA= (M2,I) *C (1,I)^A (3,1) *C(2,I) ) /Cl 
A (3,1) = (A(3,I)*C(1,I)-A(2,I) *C(2,I) )/Cl 
A(2,I) =TA 
next I 

! L is the distance from the sample to the center of 
the farthest mic, 

! S is the spacing between the microphone cenners 
! D1 is the frequency spacing of the 5420 samples 
S=.0459 

L=. 1 142 a D1=6 .25 

! calculate transfer function, impedance, reflection 
coefficient and phase change 
for 1=1 to 256 
F= (1-1) 

K=2*PI*F/C 

RH=A(2,I)/A (1,1) a IH=A (3,1) /A (1 ,1) 

AA = SIN (K*L) a BB=SIN(K*(L-S) ) 3 CC=COS (K*L) 3 DD=C0S 

(K* (L-S) ) 

IN1= BB-AA*RH 3 RN1=AA*IH 
RD1=RH*CC-DD a ID1=CC*IH 
H2 = RD1=«'RC1 + ID1*ID1 

R(I) ,RZ= (RN1+RD1-HN I^IDI) /a2 3 M (I) , I Z= (IN 1 * RD 1 -R N1 * 

ID 1) /M2 
BD= (RZ+1 ) /2 

ID=IZ/2 3 M=RD=<'RD+ID*ID 
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1120 



RPHO = -(RD/M-1) $ IRH3=-(ID/M> 

A1 (I) =SQR (RRHO=«=RHKO+IRHO*IRHO) a) ?HI(I)=-ArS2 (I R HC , RRH01 ' 
nsxt I 



! plot results 

plotter is 705 2 lorg 5 2 F0=500 2 F1=1600 



Y1=-10 2 Y2=10 2 S1=1 

disp "resistance" 2 1=1 2 L$=" 

IMPEDANCE" 

gosub 1070 2 if Z1 then 1020 
Y 1=-10 2 Y 2=10 2 S 1=2 
disp "reactance 2 1=2 2 L$=" 
IMPEDANCE" 

gosub 1070 2 if Z1 then 1020 
Y1=0 2 Y2-1.4 2 S1=.2 
disp "reflection" 2 1=3 2 LS=" 
COEFFICIENT" 

gosub 1070 2 if Z1 then 1020 

Y1=-1 30 2 Y2 = 1 80 2 SI ='45 

disp "phase angle" 2 1=4 2 L$=" 



RESISTIVE" 2 BS=" 



REACTIVE" 2 BS=" 



REFLECTION" 2 B$=" 



PHASE" 2 B$=" 



ANGLE (DEG)" 

gosub 1070 2 if not Z1 then goto 430 
disp "lowest Y" ;2 input Y1 
disp "highest Y" ;2 input Y2 
disp "increments of" ;2 input Si 
on I goto 910,940,970,1000 
! plotting subroutines 

disp "change paper and press end line";2input AS 
Y0= (Y2-Y1) /7/25 2 Y3= (Y2- Y 1 ) /7 * 1 . 9 2 Y4= (Y2- Y 1 0/7*1 .6 
2 Y5 = Y1-Y3 2 Y6=Y2+Y4 

lorg 5 2 scale - (F 1/5 ) , 1 . 16 *F1 , Y5 , Y6 2 xaxis Y1,F0,0, 

FI 2 yaxis 0,S1,Y1,Y2 

Y7=Y1-Y3/6 2 Y8=Y1-Y3/2 2 Idir 0 



for X=0 to FI step FO 2 move X,Y7 2 label VAL$ 



(X/1000) 2 next X 

move .4*F1,Y8 2 label "FREQUENCY-KHZ" 
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1220 
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1240 



for Y = Y1 to Y2 £t9p 3 1 3 move - (FI/20), Y a) label V.-.L 
(Y) 5) next Y 

ilir ?I/2 2 Y9 = (Y2-Y1 ) /3 + Y1 S EOV 2 -(FI/10), Y9 3 lab 

L3,B$ a) F9=0 

for L=1 to 256 dtep 3 

on I goto 1170,1130,1190,1200 
F2=R (L) 5> goto 1210 

F2=M(L) 5) goto 1210 

F2=A1(L) d goto 1210 
F2 = PHI (L)/PI*180 

move F9,F2 a> F9 = F9+3*D1 d imove - (FI/500) , YO/2 
idraw FI/250, 0 3 idraw 0,-Y0 3 idraw -(FI/250),) 3 
idraw 0,YO 
next L 

disp "1 for redraw, 0 to continue" input Z1 S retu 
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